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Abstract 

We investigate the late-time behavior of particle creation from an extremal 
Reissner-Nordstrom (RN) black hole formed by gravitational collapse. We cal- 
culate explicitly the particle flux associated with a massless scalar field at late 
times after the collapse. Our result shows that the expected number of par- 
ticles in any wave packet spontaneously created from the "in" vacuum state 
approaches zero faster than any inverse power of time. This result confirms 
the traditional belief that extremal black holes do not emit particles. We also 
calculate the expectation value of the stress energy tensor in a 1+1 RN black 
hole and show that it also drops to zero at late times. Some comments on 
previous work by other authors are provided. 
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1 Introduction 



As a quantum effect, particle production by black holes (Hawking radiation) was 
widely studied since the 1970s pQ 00 • It is well known that black holes emit particles 
with the same spectrum as a blackbody with a temperature T = hn/2irk, where k 
is the surface gravity of the black hole and k is the Boltzmann's constant. The 
Hawking radiation can be derived for a spacetime appropriate to a collapsing body. 
At a late stage, the collapse settles down to a stationary black hole. Since spacetime 
is asymptotically flat at both past infinity J?~ and future infinity J^ + , we can define 
the "in" and "out" vacuum states, respectively. If the two vacuum states are distinct 
from one another, particles will be detected at future infinity when the initial state is 
in "in" vacuum. If such creation takes place at a steady rate at late times, it indicates 
that those particles are produced from the stationary black hole instead of from the 
collapse phase. 

Although the standard derivations of Hawking radiation only deal with nonex- 
tremal black holes, it is generally accepted that extremal black holes have zero tem- 
perature and consequently no particles are created. However, Liberati, et al. j3] 
pointed out that the generalization to extremal black holes from nonextremal black 
holes is not trivial. One important difference is that the Kruskal transformation 
for non-extremal black holes, which plays a crucial role in computing the particle 
creation, breaks down for extremal black holes. 

The arguments in [3] are reviewed briefly here as follows. Start with the usual 
form of the Reissner-Nordstrom (RN) geometry with parameters Q and M 



where dfl 2 is the metric on the unit sphere. The tortoise coordinate r*(Q, M) is given 



ds 2 




(1) 



by 




(2) 



In the nonextremal case \Q\ < M, 




(3) 



where r± = M ± y/ M 2 — Q 2 . Define the retarded time u and advanced time v as 



U = t—T; 



V 



t + TV 



(4) 
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The well-known Kruskal transformation for the nonextremal case is 



U = -e~ RU = --\n(-U), (5) 

K 

V = e KV ^ v = -ln(V), (6) 

where U and V are regular across the past and future horizons of the extended 
spacetime. In the extremal case \Q\ = M, the right-hand side of Eq. appears to 
yield the indeterminate form 0/0. This can be fixed by setting \Q\ = M in Eq. (J2J) 
before integrating. Thus, 

n(M, M) = r + 2M (ln(r - M) - 2(r f M) ) ■ (7) 

Since k = when \Q\ — M, the Kruskal transformation (JHJ) and (jHJ) breaks down for 
the extremal case. A generalization of the Kruskal transformation to the extremal 
RN black hole is 0] 

u = -4M (ln(-U) + , (8) 

v = AM^HV)-^. (9) 

We shall show, in section 2.2, that Eq. (jHJ) defines a smooth extension. However, 
Liberati, et al. 0] actually used a simplified extension 

2M 2 , x 

« = — jp, (10) 

which, as we will show later, is not a smooth extension (Eq. (|10p is essentially the 
same extension introduced by Lake [5]). By using this extension, Liberati, et al. 
jl] calculated the Bogoliubov coefficients associated with plane wave solutions of a 
massless scalar field in a two-dimensional Minkowski spacetime with a moving mirror 
(serving as a timelike boundary) which is physically equivalent to a (l+l)-dimensional 
model of an extremal RN spacetime formed from a collapsing star. The result shows 
that the Bogoliubov coefficients are nonzero, indicating that particles are created in 
the late stages of collapse. Further calculations in [3] also show that the expectation 
value of the stress-energy-momentum tensor is zero and its variance vanishes as a 
power law at late times. The authors thereby claim that the extremal black hole does 
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not behave as a thermal object and cannot be regarded as the thermodynamic limit 
of a nonextremal black hole. 

However, the major deficiency in the analysis of [I] is the use of unnormalized 
plane-wave solutions. These kinds of solutions have been used for nonextremal cases 
[Hj [7j [H] • The Bogoliubov coefficients f3 uu i in jlj have the form 

[ V^+^cfc. (li) 

The integrand is oscillated with constant amplitude. So the integral is not well- 
defined. The result 

aujuj 1 ) (12) 

given in |3] , which was originally calculated by Davis and Fulling J5] , was obtained by 
Wick rotation, i.e., integrating along the imaginary axis. But this Wick rotation is 
unjustified since the integrand does not fall off at large radius on the complex plane. 
Since the spectrum of particle number created from the vacuum is 

<N U >= / (13) 
Jo 

and K\(z) ~ 1/z for z — > the number of particle is divergent. The authors 
interpret the infinity as an accumulation after an infinite time. The Kruskal extension 
(fTUjl was used in deriving Eq. (fTTjl. If we use the smooth extension (jHJ) instead, the 
Bogoliubov coefficients would be 

^ „ y^y 00 e -^+- e -^Hv) dv (14 ) 

and by using the same Wick rotation (also unjustified), it follows that 

~ uj'^K^m^V^JJ). (15) 

For small z } K x _± iMul {z) w |r(l - AiMoo)(^)- 1+4iMuj . Therefore, the number of 
particle in Eq.([13p is still infinite. 

To clarify this issue, our main calculation focuses on wave-packet solutions with 
unit Klein-Gordon norm. The wave packets P n euj im we will construct are made up of 
frequencies within e of Uq. They are peaked around the retarded time u = 2im and 
have a time spread ~ 2n/e. The created particle number, N ne (uo), associated with 
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the wave packet has a direct physical interpretation: N ne (uj ) is proportional to the 
counts of a particle detector sensitive only to frequencies within e of u and angular 
dependence Y\ m which is turned on for a time interval lnje at time u = Inn. Our 
calculation shows that for fixed c^o, e, I and m, N ne (u ) drops off to zero faster than 
any inverse power of n. Therefore, the traditional belief that extremal black holes 
do not emit particles is confirmed. Furthermore, if we sum N ne (uj Q ) over the integers 
n, we still get a finite result. This indicates that, even after an infinite time, the 
accumulation of particles for a certain frequency is still finite. This contradicts the 
infinite result in 0j. Our calculation is independent of choice of a specific type of wave 
packet provided that its Fourier transform is a C°° function with compact support 
on purely positive frequencies. As explained in section 2.4, we also conjecture that 
our result is independent of the details of the collapse. 

Note that two independent errors were made in First, the nonsmooth Kruskal 
extension (fTU|) was used rather than the smooth Kruskal extension (|SJ). If Eq. (|10j) were 
used in the wave-packet method, the Wick rotation used in calculating the negative 
frequency part of the wave packet at the past infinity would not be justified(See foot- 
note 1). Second, unnormalizable plane waves were used rather than normalized wave 
packets. Even if the Kruskal extension (jHJ) had been used, the use of un-normalizable 
plane waves would have resulted in the prediction of an infinite number of particles. 

We also calculate the expectation value of stress-energy tensor < T uu > related to 
the extension (|32p and find that it drops to zero as This conclusion is proved to 
be independent of the details of the collapse. From the particle flux in a wave packet, 
we find that N ne (ui )(ui )uj duj drops as fast as or faster than 1/u, which is not 
in contradiction with the \ decay rate. 

Our calculation follows the similar steps to pp. We focus on a massless scalar 
field on an extremal RN black-hole spacetime which is formed from a collapsing star. 
We start by constructing a positive frequency (relative to retarded time u) wave 
packet P neuJ0 i m at future infinity . By using the geometrical optics approximation, 
we propagate the solution back to the past infinity J 2 ^. The particle number in this 
mode can be obtained by computing the Klein-Gordon norm of the negative frequency 
(relative to advanced time v) part of P neuJ0 i m at J?~ . 
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2 Calculation of particle creation 

2.1 Construction of the wave packets at future infinity 

Our purpose in this subsection is to construct positive frequency wave packets at 
future infinity We start with the massless Klein-Gordon equation 0<p = 0. In 
the region outside the collapsing matter, the spacetime is described by the extremal 
RN metric (JTJ). Write (ft = r _1 F(r, t)Yi m (6, <fi), where Yi m (9, 4>) is a spherical harmonic. 
Then outside the collapsing matter, □</> = yields 

d 2 F d 2 F 

where 

V(r) = - (M ~ (2M 2 - 2Mr + 1(1 + l)r 2 ). (17) 



Furthermore, assume F(r,t) = g{r^)e tuJoU , where u = t — r*. Then Eq. ()16|) becomes 

2iu; </(r*) - </'(r*) + ^(r.) V(r) = 0. (18) 
As r — ► oo, r*/r — > 1 and V(r) — > 0. So <?(r*) approaches a constant and 

P.*lm = -e^ tt 15m(M) (19) 

r 

is a solution near J^ + . A wave packet with frequencies around uq and centered on 
retarded time u = 2mi can be constructed by superposing the above spherical waves 
as 

Pneuolm = ~ Zu n(u)Y lm (9 , 0) , (20) 



W(«) = ^4= r / ( ) e-^e^dw (21 ) 



where 



V 27T 

where A is a normalization constant and f(x) is a real C°° function with compact 
support in x G [—1, 1]. To guarantee P neLUO i m has positive frequencies near u > 0, we 
require < e <C cu . Let i£; ^ a = Q and it — 2nix = u. Then we may rewrite Eq. (}2*T]) 

as 

Zuim {u) = Aee**> a f(eu), (22) 
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where 



1 



f(u)e iG)X dti. 



(23) 



/2tt J- c _ 

The normalization constant A is determined by the Klein-Gordon inner product 



9J 



Taking S to be J + , Eq. (j2U) becomes 

^«9P „<9P 



;p p)^ = z J(Pv a p - pv„j>w = i 



— 2 



P^— - P^— ) r 2 dVtdu = 1 
ou ou 



Substituting flU into gives 



[^(u)z'(u) — zz'iu)] du = 1, 



(24) 



(25) 



(26) 



where we have omitted the subscripts of z^ n . Straightforward calculation from Eq. 
(|22| gives 

z(u)z'(u) — zz'(u) 

/m' 2 



= i2cj ^ 2 e 2 

+ — i2A 2 e 3 Pe 
2vr 

The solution for A is found to be 

.4 



f{Cj)e ie ™dCj 



f(u)u>e- ifS *du> 



y/ (3cu + 7e 2 ' 
where (3 and 7 are integral constants defined by 



(27) 



(28) 







c/.r 



-00 
00 



fix) 



1 
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\f(Co)\ 2 dCo 

—00 

~ />oo 

_J —00 
^|/(^)| 2 ^ 



(29) 



2— Re 
2vr 



e 4 ^/(cu)c2cu 



e- ix *f(u)Qdu 



(30) 



Therefore, Eq. (|22j) becomes 



z{u) = v e^f(eu). (31) 
2.2 Kruskal coordinates 

Kruskal coordinates will play an important role in our following calculation. Specif- 
ically, we seek a coordinate U which is a smooth function of u outside of the black 
horizon and covers a neighborhood of the horizon with U = on the horizon such 
that the metric in the coordinates (U, v, 8, <fi) is smooth on the horizon. Note that r 
is a smooth function (it is easy to check that r is an affine parameter of an incom- 
ing null geodesic). Thus, along an ingoing null geodesic with constant v, an affine 
parameter U can be taken as U = — (r — M) and U(u) is obtained from Eq. (0) by 
using r* = ^{v — u). Therefore, we have 

u = 2U - AM (ln(-U) + |Q + d, (32) 

where d is a constant. Without loss of generality, we may assume d = 0. The coor- 
dinate U defined along the ingoing null geodesic can be "carried" away by outgoing 
null geodesies. For a smooth two-dimensional spacetime defined by (U, v), there exist 
smooth coordinates (U, V) such that the metric takes the form 

ds 2 = -n 2 (U,V)dUdV, (33) 

where U = U(U). Since U must be a smooth function of U with nonzero first 
derivative along the null ingoing geodesic where U is an affine parameter, U = U(U) 
defines a smooth coordinate transformation. Therefore, we have constructed a smooth 
Kruskal extension U. Next, we wish to show that Eq. (JHJ) also defines a smooth 
extension. To distinguish, we rewrite U in Eq. (JSJ) as U'. So we have, from Eqs. (JHJ) 
and ((22) 

- 4M (in(-U') + ^j=2U- AM (\n(-U) + |0 (34) 

Since U' is obviously smooth outside the black hole, we only need to show that U' is 
a smooth coordinate in a neighborhood of the horizon, i.e., U' is a smooth function 



8 



of U and sjfr ^ around U = 0. Let £/' = [1 + Uh(U)} and substitute into Eq. JJS}. 
Differentiating both sides of (J33|) with respect to [/, we obtain 

d/i (l + (M + f/)fr) 2 
<fl7 ~ M(M -2U- 2U 2 h) 



(35) 



It is easy to check that the right-hand side of Eq. is a smooth function of (U, h) 
at U = h = 0. Therefore, according to the theory of ordinary differential equations, 
there exists a unique smooth solution to h around U = with h(0) = 0. It is also 
easy to check that ^ 7^ around [7 = 0. Therefore, [/' is a smooth extension, i.e., 
the U defined by (jHJ is a smooth coordinate. Now we show that U defined by (fTUj) 
is not a smooth coordinate on the horizon. Rewrite U in (fTUj) as U' . Then, from (|8|) 
and (fTUj) . we express C/' as 

_ Mg 

^ " M + 2U\n(-U) W 

It is straightforward to show that is divergent at U = 0. Therefore, the extension 
defined by Eq. (jl(J|) is not smooth. 



2.3 Geometrical optics approximation 

To calculate the particle creation rates at late times, we need to propagate the wave 
packet ()2()j) backward from J? + to . For simplicity, we first investigate the prop- 
agation of solution (jl9|) ; later, the propagation of Eq. ()20j) can be easily obtained by 
superposition. A part of the wave (fE?|) will be scattered by the static Schwarzschild 
field outside the collapsing body and will end up on with the same frequency [J 
and will not contribute to particle creation. We are interested in the remaining part 
which will propagate through the center of the collapsing star, eventually emerging 
to . Consider the solution (JT§j) propagating to a point x, which is very near the 
future event horizon Ti + and outside the collapsing body (see Fig. 1). The solution 
near x has a form similar to its form on 

Puoimix) ~ t(uj )e^ u Y lm (6, 0), (37) 

where t(uo) is the transmission amplitude describing the fraction of the wave that 
enters the collapsing body. Note that near the horizon, the effective frequency will 
be arbitrarily large 0. So the amplitude of Eq. (|37)l changes much slower than the 
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Figure 1: The Penrose diagram of a spherically symmetric collapsing body producing 
an extremal RN black hole ^0]- (The body is shown collapsing to a singularity at 
r = 0, but it might instead re-expand into the new asymptotically flat region; see 
[10J for further discussion of the behavior of charged shells. The behavior of the body 
after it crosses the event horizon is, of course, not relevant for our analysis.) 

phase. Consequently, the geometrical optics approximation becomes valid for the 
propagation from x back to J?~ . So the wave takes the form 

P^im = 9(t,r)e tuJoS Y lm (0A), (38) 

where S is called the phase of the wave. Each surface of constant S is a null hypersur- 
face |H] and consequently k a = V a S is the tangent to the null geodesies propagating 
in the radial direction. If we follow a light ray backward in time, it will pass through 
the center of the star and propagate to J?~ . We fix (6, <j)) and consider the family of 
radial null geodesies 7 X (A) such that for each x, 7x(A) represents a null geodesic with 
parameter A sent from radially to the collapsing star. So all geodesies in 7 X (A) 
have the same path in space but they pass through the center of the star at different 
times. The limiting null geodesic in this family lies on the future horizon. Set x = 
for this geodesic and denote it by /i, i.e., fi = 7o(A). Let x be an event lying just 
outside the horizon (see Fig. 1). According to geometrical optics [§], 5" = u(x) along a 
null geodesic. To find out the explicit form of S, let v be the Killing/afline parameter 
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coordinate at past null infinity and v = correspond to the light ray on the horizon. 
For a fixed collapse, define v on the spacetime by propagation from past null infinity 
along radial null geodesies. Let U be a smooth Kruskal extension such that it is a 
constant along each outgoing null geodesic and U = on the horizon. Therefore, a 
function v(U) can be constructed from those radial null geodesies with v(0) = 0. The 
exact form of v(U) should be solved from the geodesic equation which depends on the 
details of the collapse. However, no matter what the details of the collapse are, the 
corresponding spacetime must be smooth. Consequently, the geodesic equation is a 
smooth equation and thereby v(U) is a smooth function for U < 0. Equivalently, each 
smooth Kruskal extension U should correspond to some smooth collapse spacetime 
for which the propagation of radial null geodesies from future infinity to past null 
infinity is given by v = U(u). Thus, for the Kruskal extension defined by Eq. (|SJ), we 
have 



There is a close analog between four- dimensional spherical collapse and two- 
dimensional Minkowski spacetime with a moving mirror. The physical relations have 
been widely discussed in previous literature, e.g., jHj, [7] and jH]. We shall only il- 
lustrate the mathematical correspondence between a spherical collapse and a mirror 
trajectory. In a two-dimensional Minkowski spacetime with double-null coordinates 
(u,v), a moving mirror servers as boundary of the spacetime. If a left-moving light 
ray with constant v is reflected by the mirror, it then becomes right-moving with 
constant u. The relation between u and v is uniquely determined by the coordinates, 
(u,v), of the reflecting point on the mirror. Thus, the mirror trajectory u = u(v) 
shows how a light ray propagates after reflection. Let the left-moving light ray cor- 
respond to an ingoing light ray in a spherical collapse and the right-moving light ray 
correspond to an outgoing one. Then we see that the mirror plays the role of the 
origin of spherical coordinates. The trajectory associated with the collapse above is 
simply Eq. (J3~9*|). Such a relation will be used later to calculate the energy flux. 

The amplitude in Eq. (|55|) near J 2 ^ can be calculated by substituting Eq. ()38|) 
into □</> = 0. After neglecting the Ug term, which is supposed to be small in the 
geometrical optics approximation, we obtain 



Note that V a k a is the expansion of the congruence of radial null geodesies, which is 
equal to [3] 




(39) 



v 



2k a V a g + gV a k 



0. 



(40) 



1 dA 



(41) 
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where A is the parameter of k a and A is the cross-sectional area element. Since all 
geodesies represented by k a point radially and the spacetime is spherically symmetric, 
A oc r 2 . Hence, Eq. pUj) gives 

= 0. (42) 

d A 

Namely, gr is a constant along each null geodesic. So gr is proportional to t(u ) in 
Eq. Finally we get the solution near the past infinity, 

^ im ~| 0> w > Q) (43) 

where S'(v) is given in Eq. (J59j) . Next, superpose solutions the same way as we did on 
(refer to ([21))) and assume that t(u) varies negligibly over the frequency interval 
2e. Then, we obtain the wave packet at J?~ , 

f t(u )lz p (v)Y lm (9,(f)), v<0 . , 

ne " olm \ 0, u > 0, 1 ' 

where 



-ic;o2n7r e ^o[-4Mln(-»;)-^-] J( e {2) (45) 



and 



2M 2 

tt = -4Mln(-t>) 2wk. 



(46) 



2.4 Calculation of particle creation 

We shall show that the particle creation rate for each mode with fixed Uq, e, I and 
m will drop off to zero at sufficiently late times. So in this subsection, we treat u 
and e as fixed and consider the limit where n is allowed to become arbitrarily large. 
The expected number of particles spontaneously created in the state represented by 
a wave packet is given by jHj 

N ne (u ) = (p-,P-) KG , (47) 
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where P~ represents the negative frequency part of the solution (|44|) . The negative 
frequency is with respect to v. Since the only v dependence in Eq. (|44|) is z p (v), after 
integrating on J^~, (}4Tj) reduces to 



r>oo 

N n€ (cu ) = \t(uj )\ 2 I \z{u')\ 2 u'du\ (48) 



where 



r-0 

z{J) = I e iw ' v z p {v)dv (49) 

oo 

7 r° 2 

' / jo/-. _i,.,„o^,^ r a j\/ri„/ /„A 2aAf ' 



e 



/„ e -^ 2n V c [-4M ln(-«/o)- /(gft)^ ( 5 Q) 



OO 



is the amplitude of the negative frequency part of z p (v). Note that in z p (v), v is always 
multiplied by an undetermined factor 1/ct. By a simple rescaling, we see immediately 
that N nt (uo) is independent of the choice of a. So without loss of generality, we set 
a = 1 from now on. The difficulty in evaluating Eq. (|5()jl is the oscillation in the 
integrand. We wish to eliminate this oscillation by a Wick rotation. We shall show 
that the integral in Eq. ()50|) can be performed along the positive imaginary axis in 
the complex v plane. To justify the rotation, we need to use the following theorem, 
which corresponds to one direction of the Paley- Wiener theorem [TT] . 

Theorem 1 Let f : 1R —>■ 1R be a C°° function with support in [—1,1]. Then, the 
Fourier transform, f((), of f is an entire analytic function of ( such that for all 
k > 0, 

1/(01 - (^hclF (51) 

for all ( G C, where is a constant which depends on k. 

We shall be interested in large k. Hence, k is assumed to be large in the rest of the 
paper. Applying this theorem to f(x) defined in (j2HJ) with x replaced by eu, we find 
immediately 



;i + \eu\y 



!/(-)! ^ TFrTT^ (52) 



Thus, in contrast to the integrals considered in jlj^, the integral appearing in Eq. (joTIj) 
is convergent. The theorem also tells us that f(eu) is an analytic function of eu. 
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Figure 2: Integration contour in the f-plane 



Thus, the integrand of Eq. (|5()jl is analytic everywhere in the second quadrant except 
at the origin of the f-plane. However, we can choose a contour which goes around the 
origin along a circle with infinitesimal radius in the first quadrant. In order to apply 
Cauchy's integral theorem, we choose the closed contour as shown in Fig. 2, where 
C\ and C*2 are two circles with small and large radii, respectively. We are going to 
show that the integration in Eq. (JHUj) over the two circles is negligible. As for the 
small circle C\, we need to show that the integrand is not divergent in the second 
quadrant near the origin. From Eqs. and ([52ft . it is easy to see that the only 

possible source causing the divergence near the origin is e' /m t e °» However, this 

2 f 

. 2etiW ojqoj _ __ _ 

term is always suppressed by e~ l « in Eq. (}45j) since ujq > e. Therefore, the 
integral over C\ approaches zero. To deal with the integration over C2, we introduce 
the following lemma. 

Lemma 1 If F(z) satisfies lirnui_ ! . 00 = in the second quadrant, then 

Jo e%cz F{ z ) = when the radius of C% approaches infinity, where c> is a constant. 

The proof of the lemma is given in the appendix. From Eq. (|HU|) and (|52jl . it is 

easy to see that the condition of the lemma is satisfied. 1 Therefore, the integration 

keeping the logarithm in the u(v) is essential to make the condition satisfied. Otherwise, z p (v) 
will not drop to zero for large So the rotation does not apply to the accelerated mirror case. 
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over C*2 approaches zero. Thus, the integral in Eq. (J5U)) can be performed along the 
positive imaginary axis in the f-plane. Let x = —iu'v. Then the integral in Eq. (|5Uj) 
corresponds to the positive real axis in the x-plane, 

/~Z tD —iuj 2niT poo o / 

z(u/) = V %= / e-*e-^ Ml ^e--^f(eu(x))dx, (53) 



-iuj'y/fiujQ + 7e J 
where 

= -4Mln( — ) H 2nvr. (54) 

ioj' x 

Using the fact that Im[u] = 2Mvr + and | e -*-o^»»(^r)| = e -2M™ when x ig 

real, together with Eq. (J53J). we have 

1 1 ;I " o;V^o + 7e io (l + |e M (x)|) fc 1 ; 

To proceed, we first derive a lower bound for \u(x)\ at large n . Start with 

(2n7r + AM H^)) 2 + ( 2M7r + ^7^") • ( 57 ) 



|m(x)| 2 



Let v = ^7 and define 

{ 2M 2 \ 2 

h(y) = \u(x)\ 2 = (2mr + AMlnyf + \ 2Mn + J (58) 

To find the minimum, we solve h'(y) = 0, which gives 

- — (8M 4 + 8M 3 vn/ - 16Mmn/ 2 - 32M V In y ) = 0. (59) 

Obviously, y — is not a solution where /i(y) achieves its minimum. When n is a 
large number, the solution to (J59j) must be at small y. So the approximate solution 
is 

y° = -n^ ( 6 °) 

Vnv 27r 
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and 



(61) 



It is easy to check, by computing the second derivative, that h(yo) is a minimum for 
large n. So we have an important inequality 



\u(x)\ > 2rar 

for large n. Then it follows immediately from Eq. (|56|) that 



\z(u')\ < 



y/eC k e 



-2A/7r(w -e) 



1 



u'y/fiuo + 7e e k (2mr) k J 



e e 



dx. 



(62) 



(63) 



However, this bound is not good enough for all u' since we must integrate \z(uj')\ 2 uj' 
over all to' and the bound in Eq. (J53j) will lead to a divergence at small uj' . To avoid 
this divergence, we need to investigate the bound in Eq. (JSTIj) more carefully. Denote 
by G{ui') the integral in Eq. (|56[). i.e. 



First, rewrite |w(a;)| 2 in Eq. fjoTj) as 



00 e~ x e 



2M 2 u'(u()-i) 



-dx. 



\eu\x) 



(64) 



|m( x )| 2 = (F n (u') + AMlnxf + 2Mvr + 



2MV^ 2 



where 



FJu/) = 2nvr -4MhW. 



(65) 



(66) 



In the following discussion, we choose ficl and consider the frequency range u/ < Q. 
For ?!>1, F n {uj') is a large positive quantity. Substitute Eq. (JfiHjl into Eq. (JH3J) 



GV) 



(i^u/) + 4Mlna;) 2 + (2Mvr + ^f^) 2 



fc/2 



(67) 



I will evaluate the bound for this integral in three domains of x. Choose a such that 
o€l and an 1. 
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(1)D1 = { e -«W)/4Af < x < e *(-')/4M and x > wjV } 
The integral in this interval is approximately 

e e 



G(u',Dl)= — dx. (68) 

Jdi e fc [F„(u/)J 

Since Dl is a subset of (0, oo) and the integrand is always positive, we have 

G(u',Dl) < r / e~ x e ^^dx 

[eF n (u/)]* Jo 

= r 2 ^ Nlfc V^ 2 (^o - ^V^i (2v / 2 V / ^ 2 (^o - e)^) (69) 
[e^n(^)] V ' 

where K\ is a modified Bessel function. For z — > 0, we have if 1(2) — ■> -j; There- 
fore, for small a;', (|69j) takes the form 

(2) ^ = {x < ^} 
The integral in this range is 







Note that aF n {u') < gives 



< e -a-fn(w ){fJ -e) ^2) 

Replacing e _:c by 1 and using Eqs. ()62|) and (J72j) to replace the corresponding terms 
in the integrand of Eq. (|7T|). we obtain the bound for Cr(u/, -D2), 

G(u', D2) < -±—e- aF ^^^^-. (73) 
v ' ; - (e2mr) fc aF n (u') V ; 

(3) D3 = {x E (0, oo)|x <£DlU D2] 

Obviously, D3 is a subset of all positive x satisfying aF n (u f ) < |4Mlnx|. Then, it 
follows from Eq. (J64j) that 

e * , f e x e 



Gfu/, .D3) < / . — ; . ., dx + / _ , ,, : — i—^—. dx. (74) 

17 



2M 2 J(u -e) 



Replace \u(x)\ by 2mr and e x e * by 1 in the first integral. Replace 

_ JM S u'(lafl-e) 

by 2nir and e * by 1 in the second integral. Then we obtain the bound for 
G(a/, D3) 

Combining (fTUj) . (J75J) and (|75j) . we have the bound for G(u/) for a/ < fi, where ficl, 



[eF„(o/)f 02n7r) fc aF n {u') 

+ t rr e + t rr e • ( 

(e2nii) k (e2n%) k 

Since F n {u') can be arbitrarily large, we only need to keep the first term on the 
right-hand-side of Eq. (|76p. This fact reveals that the integration in Eq. (joT|) is 
approximated by taking away the x-dependent terms in the denominator of the inte- 
grand. Thus, the bound for |£(u/)| in ()56|) at small uj' is 



%'<Q)1<^= : = r _ (77) 



^C k e- 2M ^°-^ 
For u/ > fi, we simply use the bound (J63J) . Then 



v^Cfee- 2 ^^-^ 2V2 



< V ,V— - ^-^VM 2 K-6KK 1 (2v / 2 v / M 2 K-eK) (78) 

Now we are ready to compute the particle creation rates at late times. It follows 
from Eqs. (JZZj) and (JZHI) that Eq. (gHJ) is bounded by 



2 cu/ fc 



(/9w + 7e)e 



2k 



-dw' 



J uj' (2mr -4MbW) 2fc 

(79) 



+ 4 r M ^ " £) ^ 2 (2v / 2 V / M 2 K - 6)o/)do/ 
n Jn u 

Evaluating the first integral gives 



AM 1 4M 1 

(2k- 1) (2nvr-4Mln^) 2fc - 1 ~ (2fc - 1) ^wr) 2 *- 1 
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The second integral in Eq. (f79*|) is convergent since K\(z) ~ \/iz e ~ Z ^ or ^ ar S e 
\z\ [I2|. Since n represents time, we conclude that the particle creation rate for any 
mode decays with time faster than any power law. Furthermore, by summing over n 
from any positive integer to infinity, the bound in Eq. (|T9~j) is still finite. This means 
that, as mentioned in the introduction, the accumulation of particles after an infinite 
time is finite. This conclusion is derived from a particular Kruskal transformation (JHJ), 
which corresponds to a particular process of collapse. A different process of collapse 
will give rise to a different Kruskal coordinate U', which is a smooth function of U. 
We conjecture that our conclusion that the particle creation rate decays with time 
faster than any power law is independent of the choice of the Kruskal extension, i.e., 
independent of the details of collapse. As evidence in this conjecture, one can check, 
following similar steps, that the smooth extension (J32J) also gives the same result. 



3 Stress energy tensor 

In two-dimensional Minkowski spacetime, the renormalized energy flux in a spacetime 
with a moving mirror boundary in the "in" vacuum state is |H| 



47T 



1 f P " 



4 \p' 



lp'" 
6 p' 



where p = v = p(u) is the trajectory of the mirror. According to the discussion in the 
last section, p is exactly a Kruskal extension U which describes a particular collapse. 
The corresponding trajectory of Kruskal extension (JHJ) is thereby 



u 



-AM ln(-p) 



2M 2 
V 



Straightforward calculation yields 

< T >- 

which approaches 



p 3 (p - 2M) 
48ttM 2 (M -2 P y 



(82) 



if 



(83) 



for large u. Therefore, < T uu > decays as 1/u 3 . 
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The quantity < T uu > can be estimated from the particle flux by adding up all 
frequency modes at certain time 



CO 



< T uu >~ / N ne (uj )u duj . (84) 



o 

This formula, as discussed in jH], is a naive energy-particle relation. It is correct 
only when particles emitted in different modes are not correlated. However, it is 
worth comparing this naive < T uu > with the one in Eq. ()83)1 . If the naive one is 
smaller, it may indicate some serious problems in our calculation of N ne (uo). From 
the estimation of N ne (u ) in the last section, Eq. suggests that the naive < T uu > 
also should decay to zero faster than any inverse power of u. So the result 1)83)1 may 
seem to be inconsistent with the particle creation results. However, in the last section, 
we treated o>o and e as fixed while allowing n to be arbitrarily large. But Eq. (1841) 
requires us to sum over all modes at a fixed time n. So we need to reevaluate the bound 
of N ne for arbitrarily small uq. In this subsection, we focus on the (l+l)-dimensional 
RN black hole formed by collapse because our purpose is to check the consistency 
of ([83)1 which is computed in a 1+1-dimensional spacetime. The calculation will 
be parallel to our four dimensional case in the previous sections. One important 
difference is that the transmission amplitude t(uo) has unit magnitude due to the 
fact that a (l+l)-spacetime is conformal to Minkowski spacetime and therefore the 
outgoing wave packet will not be scattered when it is propagated backward in time. 
So the bound in Eq. (f79|) still holds for the (l+l)-dimensional case except | t(a^ ) I = 1- 
However, for our present purposes, the bound in Eq. 1)56)1 becomes inappropriate since 
|ew(x)| ^> 1 is not always true for arbitrarily small e. So we stick to Eq. 1)5*5")) and 
follow similar arguments. Denote by H(u') the integral in (|64j) . i.e., 



„ 2M 2 u'(u n -e) 

00 e~ x e 



Jo (i + Mx)|)* 



Consequently, the corresponding changes in Eq. 1)76)1 become 



H{J) < - k + - L^ e -^K)(-o-e)^!^ (86) 

V ^ " [l + eF n (u')] k (l + e2nvr) fc aF n (u') V ; 

1 aFnK^ ) 1 _„ — AU 

+ e 4M H e 

(1 + e2nir) k (1 + e2nir) k 

The last two terms in the bound are still negligible due to the exponentials. Unlike 
before, we are not certain whether the second term is much smaller than the first 
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term since uq — e can be arbitrarily small now. So we keep both of the terms. In 
order to perform integrals easily, we replace the exponential and u' in the numerator 
by 1 in the second term. Thus, 

1 1 2M 2 

[U) ~ [I + eF n {u')f + (1 + e2nnr aF n (u>) (8?) 

Then (J55J) becomes 

lv , mi ^ v^C fc e- 2M -^~ s ) ( I 1 2M 2 

l^<^)l < ,., 7zf= =— 77- 777 771 + 



;V^ + 7e \ [1 + eF n (u/)] fe + e2mr ) fc «^(^') 
The modification for z(u>' > Q) in Eq. (|78jl is 



< 



|z(a/ > Q)| 



w'Vjfwo + 7e (! + e2mr) fe 
Thus, the bound on the particle number from is 

A^„e(^o) 

< 2eC 2 e- 4M ^°- e ) /" n 1 / 1 V 1 / 1 2M 2 V d 

~ (3u + 7 e i Z7 yi + eF n (a; , )]V + ^ V (1 + e2mr) fc aF n (o/) , 

e ^2g-4M7r(a;o-e) c /-co ^2 

+ 



/3o;o + 7e (1 + e2nvr) 



2A: 



f° M (tJ0 - e)J K\ (2T/2y/M*{uo-e)u>) du>, 
Jn v ^ ' 

(90) 



where we have used the inequality (A + B) 2 < 2(A 2 + i? 2 ). The integration over (0, Q) 
gives 

2eC|e- 4M "^°- e ) ( 1 1 M 3 



+ 777777777-7777777^), (91) 



(3ujq + 7 e V ( 2fc - l)4Me (1 + 2n7re) 2fc ~ 1 2mra 2 (l + 2n7re) 



where we have neglected the ^-dependent term since it is small compared with n. To 
estimate the integration over (O, oo), we first change the integration variable from u' 

to x = (lu - e V'. Thus the integral becomes M 2 J^ o _ e) K\ (2\f2*JWx\ dx. We 

then split the integral into two terms as 



M 2 [ K\ (2V2VM 2 x) dx + M 2 [ K\ (2v / 2v / M 2 x) dx, 



(92) 
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where b > — e) and bM 2 <C 1. Thus, we can use the approximate form of 

Ki(z) ~ 1/z again for the first term. The second term is just a constant and negligible 
compared to the first term when uq is taken small enough. Therefore, the integral is 
approximately ^7=- ln[f2(u; — e)]. Together with a coefficient, the second integral in 
Eq. (JIID yields 

eC 2 ~2M^ g _ M 2 

— t t — i=ln\Q(u - e) . (93) 

As discussed above, the main contribution to < T uu > comes from the low-frequency 
integration J Q N(u>o)uJoduJo, where 5 is a small constant. Since it is assumed e < u , 
we can no longer treat e as constant when performing the integral. For convenience, 
we choose e to be a small fraction of u . It follows from Eqs. (}9*Tj) and (J9llj) that 

f s K 

/ N ne (u )u duj < — , (94) 
Jo n 

where K is a constant. Thus, although the particle creation rate in each individual 
mode goes to zero faster than any inverse power of time, the energy flux estimated 
from the particle creation rate decays only as 1/u on account of contributions from 
the very low frequency modes. The bound in Eq. (pM|) decays more slowly with time 
than Eq. (|%3*j). This difference may be due to the fact that the bound (JiMj) is not 
sharp enough. An alternative possible explanation is that, as pointed out by Davies 
and Fulling [B] , the relationship between particle fluxes and energy fluxes can be more 
subtle than as a result of destructive interference between different modes. At 
this stage, we do not know if there exists interference since we only have the upper 
bound on the number of particles. 

One can also ask if the 1/u 3 decay rate of energy flux is independent of the details 
of the collapse. Since each collapse corresponds to a smooth extension, let us consider 
another smooth extension q and the function q = q(u), which gives the related energy 
flux, 



1 

Air 



1 (q"V lq'" 

4\7/ ~qV 



(95) 



Let q = g{p). From q{u) = g{p{u)), we derive the relationship between < T uu > q and 
< T uu > p (the energy flux associated with p) , 



< T uu > q =< T uu > p +[p {u)\ A87Tg >( p y • 1 > ()) 
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Since p and q are two smooth coordinates, the derivatives of g(p) must be finite and 
g'{p) 7^ 0. From Eq. (|81|). [p'{u)} 2 goes as l/w 4 at late times. Therefore, the second 
term on the right-hand side of ()96|) is negligible compared to the first term which 
goes as 1/u 3 . Therefore, we showed that the 1/u 3 decay rate is invariant under a 
smooth change of mirror trajectory. This fact agrees with the fact that the late-time 
radiation is independent of the details of collapse. If one applies the extension (jTUjl 
to compute the energy flux, as done in the energy flux would be identical zero 
at late times. This is because the extension (|10|) fails to be a smooth one (see the 
comment at the end of section 2.2) at late times (around {7 = 0). 

4 Conclusions 

We have calculated the number of particles created from a massless scalar field on an 
extremal RN black hole spacetime formed by collapse. We found that for each mode 
associated with a wave packet, the rate of particle creation drops off to zero faster 
than any inverse power of time at late times. Consequently, even after an infinite 
time, the number of particles detected by a detector sensitive to a certain frequency 
is finite. This result confirms that extremal black holes do not create particles. In 
the (l+l)-dimensional case, the stress-energy flux falls off as 1/u 3 . This result is not 
in contradiction with the much more rapid decay of each mode, since the very low 
frequency modes do not achieve their asymptotic decay rate until very late times. 
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Appendix: Proof of Lemma 1 

Let R denote the radius of C 2 and Ir = J c e %cz F{z). Substitute z by Re w . Then 



Acknowledgments 




(97) 
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Since limui^oo 1^(^)1 — > 0, for any e' > 0, we can find R such that 



\Ir\ < e'R I e~ cRsine de 

h 



< e'R 2 e- cRcose d9. (98) 
Jo 



In the range [0, |], 



Therefore, 



cos9 > 1 - -9 

71 



< e'R / 
Jo 



e -cR(i-ie) d9 



Therefore, 



~ 2c v 1 



lim 1/^1 = 0. (99) 

R-^oo 
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